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The problem of the evaluation of the two-photon decay width of excited states in hydrogen is considered. Two 
different approaches to the evaluation of the width including cascades channels are employed: the summation of 
, the transition probabilities for various decay channels and the evaluation of the imaginary part of the Lamb shift. 

' As application the decay channels for the 3s level of the hydrogen atom are evaluated, including the cascade 

' transition probability 3s — 2p — Is as well as "pure" two-photon decay probability and the interference between 

both channels. An important role should be assigned to the "pure" two-photon probability in astrophysics 
■ context, since processes of this kind provide a possibility for the decoupling of radiation and matter in the early 

' Universe. We demonstrate the ambiguity of separation of the "pure" two-photon contribution and criticize the 

, existing attempts for such a separation. 
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I. INTRODUCTION 



During recent years the two-photon decay processes in hydrogen have attracted special attention due to new and very accurate 
observations of the cosmic microwave background temperature and polarization anisotropics IHSl • In view of these observations 
O it becomes important to understand the hydrogen recombination history with high precision. In the early Universe the strong 
Lyman-alpha 2p— Is transition did not permit the atoms to remain in their ground states: each photon released in such a transition 
^ in one atom was immediately absorbed by another one. However, due to the very weak 2s — Is two-photon decay process the 
Q radiation could finally decouple from the interaction with matter and thus permit a successful recombination. The role of the 
" ^ 2s — Is two-photon decay was first established in QHl- The other two-photon channels, i.e. ns — Is, nd — Is decays were also 
discussed in H-llst]. At the present level of accuracy reached in the astrophysical observations these contributions also become 
f~i important. 

Qh' a crucial difference between the decay of the ns (with n > 2, similar for nd) and the 2s levels consists in the presence of 
the cascade transitions as the dominant decay channels. In case of the 2s level the cascade transitions are absent. Since the 
cascade photons can be effectively reabsorbed, the problem of separating the "pure" two-photon contribution from the cascade 
^ contribution arises. An interference between the two decay channels, i.e. "pure" two-photon decay and cascade, should also be 
\^ , taken into account. 

• A similar problem did arise within the theory of the two-electron highly charged ions (HCI) ioll- lfTlll . Drake evaluated for 
the first time the two-photon ElMl transition in the presence of a cascade transition in He-like uranium ion (Z = 92). Later 
Savukov and Johnson preformed similar calculation for a variety of He-like HCI (50 < Z < 92) [lOJ. In l^d^ the "pure" two- 
photon contribution was obtained by subtracting off a Lorentzian fit for the cascade contribution from the total two-photon decay 
frequency distribution. The existence of the interference terms was recognized in |9, 10], but only approximately included in the 
Lorentzian fit as an asymmetric deviation from the pure Lorentzian. A rigorous QED approach for the evaluation of two-photon 
decay width with cascades was developed in ifTlll (see also |0)- This approach was based on the standard evaluation of the 
decay probability as a transition probability to the lower levels. In case of cascades the integral over emitted photon frequency 
distribution becomes divergent due to the singular terms, corresponding to the cascade resonances. To avoid such a singularity, 
the resummation of an infinite series of the electron self-energy insertions was performed in liTlll . This resummation converts 
into a geometric progression and in this way the self-energy insertion (and the level width as its imaginary part) enters in the 
energy denominator and shifts the pole on the real axis into the complex energy plane, thus making the integral finite. With this 
approach F. Low first derived the Lorentz profile from QED 1 1 30 . In 13, JO] the level widths in the singular energy denominators 
were also introduced, though without special justification. Similarly, i.e. by introducion of the level widths in the singular energy 
denominators the two-photon decay of ns and nd excited levels in hydrogen was evaluated most recently in the astrophysical 
papers fllB. 

In []JJ the ambiguity of the separation of the "pure" two-photon decay and cascades was first revealed for HCI, where it 
was shown that the interference terms can essentially contribute to the total decay probability. The full QED treatment of the 
two-photon decay process with inclusion of cascades and interference terms was performed in fT?, T?] for the hydrogen atom. 
However, in |14, 15] the ambiguity of separation of the "pure" two-photon decay was neither emphasized nor demonstrated 
explicitly. This will be subject of the present paper. Our numerical results for 3s-level in hydrogen are in agreement with the 
most accurate recent calculations in 01411 (see Section 1 1 below). However, our results disagree strongly with the value obtained 
in lfT6ll . The result in lfT6ll follows from the "alternative" approach to the evaluation of the two-photon decay, developed by 
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Jentschura iflTi [Tsi [l9|] . This "alternative" approach is based on the evaluation of the imaginary part of the two-loop Lanb shift 
for 3s-level. The idea to use the imaginary part of the Lamb shift for calculating the radiative corrections to the one-photon 
transitions is due to Barbieri and Sucher lEoll and is definitely adequate for this propose. Later this approach was used in ll2lll 
for HCI. Still, the extension of this approach to the evaluation of the two-photon decay width, when the cascade transitions are 
also present, is more involved and requires special care. In L17> .18, 191 it was claimed that the singularity in the integration over 
the photon frequency distribution for the two-photon decay is absent and thus the finite integral represents directly the "pure" 
two-photon decay rate. 

Also this statement forced us to perform a careful revision of the evaluation of the two-photon decay width via the imaginary 
part of the second order Lamb shift, which will be presented below in Sections 6-10. In our approach we employed the con- 
sequences of the "optical theorem" for the ^-matrix and the Gell-Mann and Low adiabatic formula |22] for the evaluation of 
the energy level shift via S'-matrix. The results of our analysis, unlike Jentschura's derivations, reveal that the evaluation of the 
two-photon decay width via the imaginary part of the Lamb shift provides identical expressions as the standard QED description 
via summation of transition probabilities. The integration over the emitted photon frequency for the two-photon decay with cas- 
cades remains divergent and requires the introduction of the level widths in the singular energy denominators (resummation of 
QED-radiative corrections). The detailed derivations employing both methods (summation of transition probabilities and evalu- 
ation of the imaginary part of the Lamb shift) are given and numerical results for the 3s-level in hydrogen are provided. These 
results demonstrate clearly the ambiguity of the separation of the "pure" two-photon and cascade contributions with accuracy, 
required for the modern astrophysical investigations (i.e. 1%). 

Our paper is organized as follows. In Section 2 we shortly review the standard derivation of the excited atomic level one-photon 
width via the summation over transition probabilities to the lower levels. This presentation of formulas, though available in text- 
books on quantum electrodynamics, are necessary for the further comparison with the results obtained by another approaches. 
It is for the readers convenience to make the paper self-contained and allowing to follow our derivations without consulting 
additioanl literature. The same is done (even in a more compressed way) in the Section III for the two-photon decay width. In 
Section 4 we repeat (again very briefly) the standard understanding of the situation for the two-photon decay with cascades and 
present the formulas employed for the appropriate calculations in the numerous works on the subject, beginning from [9]. For 
justification of this procedure we follow 1 11 ] . In Section 5 we evaluate the one-photon decay width via the imaginary part of the 
Lamb shift. This is again a well-known result but it will be also necessary to confirm our further derivations. In this Section we 
follow derivations of fe l in a more condensed way. 

In Section VI we evaluate the imaginary part of the energy level shift employing the adiabatic formula of Gell-Mann and Low. 
The key idea of our approach - the application of the "optical theorem" - is formulated in Section VII. Actually, we employ 
here the unitarity of the adiabatic S'-matrix. The formulated approach is applied to the evaluation of the one-photon width in 
the Section VIII, to the evaluation of the two-photon width (without cascades) in Section IX and to the evaluation of the two- 
photon width (with cascades) in Section X. In these sections we compare our derivations with the standard QED derivations and 
analyse the discrepance with Jentschura's results in Section X. In Section XI the numerical results for the 3s-level in hydrogen 
are given and the ambiguity of separation the "pure" two-photon and cascade contributions is demonstrated. Finally, in Section 
XII we formulate our version of the problem of the cascade separation in astrophysics and present our recommendations for the 
resolution of this problem. Relativistic units h ^ c = 1 are employed. 



II. ONE-PHOTON DECAY WIDTH VIA SUMMATION OF TRANSITION PROBABILITIES 

A. Transition probability 

The first-oder matrix element of the 5-matrix describing the one-photon emission in a one-electron atom (see Fig. 1) is given 

by 

{A'\S^^^\A) = e J d^x^PA'{xh^A;{x)^PA{x) . (1) 

Here S^^^ is the first-order S'-matrix, e is the electron charge, ipAix) — tpA{r)e^^^^*', tpAir) is the solution of the Dirac equation 
for the atomic electron, Ea is the Dirac energy, ^a' — '4'a'^o is the Dirac conjugated wave function with tp\, being its Hermitian 
conjugate and 7^ = (70, 7) are the Dirac matrices. The photon field Af^{x) appears in terms of eigenmodes of the form 

4"^'^^^) = y^e^e^t^— *) = y^e(^)e--4-^^)(f), (2) 

where e^^^ is the photon polarization 4-vector, k = (k, uj) is the photon momentum 4-vector (k is the wave vector, = |fc| is 
the photon frequency), x = (r, t) is the coordinate 4-vector (r, t are the space- and time-coordinates). 
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Describing emitted (real) photons in the Coulomb gauge implies the transversality condition 



(3) 



where e '^'^^ is the 3-vector of the photon polarization. Then integrating over the time variable yields 

V <^^'^ / A'A 



27re 



A'A 



Siu-EA+EA') . 



(4) 



Here a — 707 are the Dirac matrices and (. . . )a'A denotes the spatial matrix element {A'\ . . . \A) with the wave functions 
ipA'i^ = {A'lr) andi/jAir) = (r| A), respectively. 
The transition amplitude Ua'A is defined as 



{A'\S^^^ \A) = -2tt i5 {uj-Ea + Ea' ) C/^V 



A'A ■ 



(5) 



since the square of a (5-function is actually 



The transition probability cannot be defined straightforeward as 

meaningless. The standard way to overcome this technical diffculty is to express one of the two (5-function, which arises in Eq. 
O after integration over times t as the Fourier integral 



00 



(6) 



by a representation 5t{E). The latter is introduced by restricting the integration to the finite time interval {—T/2, +T/2) with 
the result 



SriE) 



1 

2^ 



T/2 



(7) 



-T/2 



Multiplying 5t{E) by the second 5-function S{E) results in the substitution 

S{E)St{E) = 6{E)5t{0) = S{E) 



T 
2^ 



Thus, the probability of the process appears to be proportional to the observation time interval T. It is natural then to introduce 
the transition probability per time unit (transition rate) by definition 



Wa'A = lim ^ 



{A'\S^^\T)\A) 



= 2ti 



rjil) 
'-^ A'A 



5{uo-Ea + Ea') . 



(8) 



We will compare this definition with another one in Section 8. 

If the final state belongs to the continuous spectrum (due to the emitted photon in our case) the differential transition probability 
should be introduced 



dWA'A{k,e) = 2tt 



U 



(1) 



A'A 



5{u}-Ea + Ea' 



dk 

(2^ 



(9) 



where dk = d'^k = d^lpduau'^ . Integration in Eq. (|9]l over w gives the probability of the photon emission with polarization e in 
the direction v = k/oj per time unit (and solid angle dfl;; = dv): 

c2 



dWA'A = -^^A'A 



((e*a)e 



-ikr 



A'A 



dv , 



(10) 



where loa'A = Ea — Ea' ■ Note, that due to the presence of (5-function formally we should integrate over the frequency interval 
{—00, 00), while physically the frequency of real photons changes within the interval (0, 00). The total transition probability 
follows from Eq. ( fTOl l after integration over angles and summation over the polarization: 



Wa'A = ^^A'aJ2 J dv\[{e*a)e-'~^^^ 

e 

The result of this summation and integration will be carried out in Section 5. 



A'A 



(11) 
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B. Nonrelativistic limit 

For the atomic electron the characteristic scales for |r| and ~ oj are: \f\ ^ 1/maZ, uj = Ea' — Ea ^ m{aZY, where m 
is the electron mass, a is the fine structure constant, Z is the charge of the nucleus. Then in the nonrelativistic case, in particular 
for the hydrogen atom (Z — 1), the exponential function in the matrix element in Eq. (fTTT l can be replaced by 1 , since kf ~ a. In 
the nonrelativistic limit the matrix element involving the Dirac matrices a (electron velocity operator v in the relativistc theory) 
can be substituted hy p/m, where p is the electron momentum operator Then Eq. ( fTTT ) takes the form 

Wa'A = ^^A'aJ2 I di^K^i^A'Al' . (12) 

e 

where the notation {■■■)a'A now also implies evaluation of the matrix element with Schrodinger wave functions. 
Performing summation over the polarization with the help of the known formula 

F^^{e*a)ieb)^{axi7){bxi7) (13) 

e 

with a, b being two arbitrary vectors, we arrive at the standard expression for the one-photon probability in the "velocity" form 
after integrating over j7: 

WA'A^^^U;A'A\{p)A'Af ■ (14) 



The so-called "length" form 
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Wa'A = 7:^A'A 



id) 



A' A 
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involving the electric dipole moment operator d = er of the electron can be obtained via the quantum mechanical relation 



(15) 



l^A'A{r)A'A = —ip)A'A ■ (16) 

m 

Eq. ( fTsT i reveals that in the non-relativistic limit the atomic radiation is essentially the radiation of the electric dipole. 
In the nonrelativistic limit the total one-photon width of the atomic level A results as 



E^,<Ea 



(17) 



where the summation runs over all levels A' with energy Ea' lower than that of the level A (provided the transitions A ^ A' 
are allowed in the nonrelativistic limit). 

III. TWO-PHOTON DECAY WIDTH VIA SUMMATION OF TRANSITION PROBABILITIES 

The two-photon transition probability A ^ A' + 2-f corresponds to the following second-order 5-matrix elements (see Fig. 

2) 

{A'\S(^^\A) =e^ J d^xid^X2 {^A'{xlh^nA;,{xl)S{x^X2h^,A*^,{x2)M^2)) , (18) 
where S'(a;ia;2) is the Feynman propagator for the atomic electron. In the Furry picture the eigenmode decomposition reads (e.g. 



(l-iO)+C^l' ^ ' 
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where the summation in Eq. ( fT9b extends over the entire Dirac spectrum of electron states n in the field of the nucleus. Using 
again Eqs. (|5]l and ^ for the two-photon transition and integrating over time and frequency variables in Eq. ( fTSl l for the sum of 
the contributions of the both Feynamn graphs (see Figs. 2a, 2b), we find 



dWA'A = 2n5 {Ea - Ea' -uj- uj') 



rji-2) 



2 dk dk' 
(27r)3 (27r)3 ' 



(20) 



^A'A — 



E 



A'n ("^^e'.fc') 

E„ -Ea+oj' 



nA 



E 



e -k' J A'n \ '^^t'/ n/ 
En — Ea + 



(21) 



with the abbreviatory notation A^j: — ee^^"^ . 

In what follows, we will be interested in the decay width of the ns levels {A = ns ^ A' = Is) in hydrogen. In this section 
we focus on the case n = 2, when the cascades are absent. In the nonrelativistic limit, after the integration over frequencies 
Lo' , over photon directions dv, dv' and summation over all polarizations e, e ' is performed, we obtain for the photon frequency 
distribution: 



dWn 



27^^ 



(22) 



n p 



En'p En 



■ LO 



(23) 



{Rn'i'\r\Rni) = J r^Rn'v{r)Rni{r)dr, 





(24) 



where = Ens — Eis, Rniir) are the radial part of the nonrelativistic hydrogen wave functions, and Eni are the hydrogen 
electron energies. Here we have used again the quantum-mechanical relation Eq. ( fTSI ): Eq. ( |22] | is written in the "length" form. 
The decay rate for the two-photon transition can be obtained by integration of Eq. ([T]) over the entire frequency interval 



Wns.ls 



wo 

I j dWns,ls{i0) 



(25) 



In case n = 2 the cascade transitions are absent, the frequency distribution Eq. (I22b is not singular and the integral Eq. ( |25l ) is 
convergent. 



IV. TWO-PHOTON DECAY WITH CASCADES 



In case of the cascade transitions {n > 2), some terms in Eq. ( |23T l become singular and the integral Eq. ( |25] ) diverges. This 
divergency has a physical origin: an emitted photon meets the resonance. So the divergency can be avoided only by introducing 
the width of this resonance. This situation was studied in jl ill for the HCI. The same recipe can also be used in case of the 
hydrogen atom. Following the prescriptions given in [TUl we separate out the resonant terms (corresponding to cascades) in the 
sum over the intermediate states Eq. ( |23] ) and apply Low's procedure 1 13] for the regularization of the corresponding expressions 
in the vicinity of the resonance frequency values. Practically this leads to the apperance of the energy level widths in the energy 
denominators. Then the Lorentz profiles arise for the resonant terms in the expression for the probability. However, the Lorentz 
profile is valid only in the vicinity of the resonance and cannot be extended too far off from the resonance frequency value. As 
for any multichannel processes such a separation is an approximate procedure due to existence of the interference terms. 

The integration over the entire frequency interval [0, cjq] in Eq. ( |25] l should be split into several subintervals, e.g. 5 in case 
of the two-photon emission profile for the 3s-level decay, see Fig. 3. The first interval (I) extends from = up to the lower 
boundary of the second interval (II). The latter one encloses the resonance frequency value oji — E^s — E2p- Within the interval 
(II) the resonant term n = 2 in Eq. ( |23] ) should be subtracted from the sum over intermediate states and replaced by the term 



with modified energy denominator This modified denominator is E2p — E; 



3s 



UJ 



^r, where T = T 



2p 



The third 



interval (III) extends from the upper boundary of interval II up to the lower boundary of the interval (IV), the latter one enclosing 
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another resonance frequency value lu2 = — Eis- Within the interval (IV) again the resonant term n = 2 in Eq. (l23T l 
should be replaced by the term with modified denominator — Eig — uj — ^2^"^^ Finally' ^ fifth interval (V) ranges from the 
upper boundary of the interval (IV) up to the maximum frequency value ujq. Note, that the frequency distribution dW-is.is iy>) is 
symmetric with respect to = a;o/2 with a 1% accuracy (the asymmetry is due to the difference between F — + and 
r2p, respectively). The discussion on the choice of the size of the intervals (II) and (IV), which defines also the size of the other 
intervals, as well as on further approximations will be postponed until Section 1 1 . 

In both Section III, IV the total width of the levels r„s is defined by the sum of the one-photon and two-photon transition 
rates to the lower levels. The cascade transitions yield the dominant contribution to r„s except for the case n — 2. 



V. ONE-PHOTON DECAY WIDTH VL4 THE IMAGINARY PART OF THE LAMB SHIFT: DIRECT EVALUATION 

In this Section we will briefly recall the well known derivation of the one-photon level width arising from pure one-photon 
transitions via the imaginary part of the Lamb shift. This may serve as a lucid introduction to the further derivations with 
employment of the adiabatic 5'-matrix theory. Let us consider a (free) one-electron ion or atom in the excited state \ A) interacting 
with the vacuum |0-y) of the quantized radiation field (with no additional external electromagnetic fields present). The initial 
state of the total system "atom h- radiation field" |/) = |yl)|0-y) = \A, 0^) is described as pure number eigenstate. Since the 
one-loop vacuum-polarization contribution to the Lamb shift is real, the pure one-photon width of the excited level A is given 
by the imaginary part of the one-loop electron self-energy contribution to the Lamb shift for the level A (see Fig. 4): 



A^;^^^ = Rel\E^l^ + ilml\E^l^ = Reb^E^l^ ~ . (26) 



Note, that the superscript at AE^'^^ refers to powers of the coupling constant e, while the superscript at r{^^ implies the number 
of emitted photons. The second-order S'-matrix element, which corresponds to the Fig. 4, looks like 

{A\S^'^^\A) = J d'^xid*X2{-ipA{xi)jf^iS{xiX2hp,2i^A{x2)) Df_,^f_,.^{xiX2) , (27) 

where D^i^j {^1^2) is the photon propagator in the Feynman gauge: 

00 

D^,^.,ixiX2) = ^^-^ [ d,^,e'-i(*i-*2)+^l-iki.^ (28) 
Zm ri2 J 

—00 

ri2 = It'i — 1- For the evaluation of this energy shift corresponding to an "irreducible" Feynman graph (i.e. such graphs cannot 
be cut into subgraphs of lower order by cutting only the electron lines; the graph Fig. 3 belongs to the "irreducible" ones) the 
following formula can be used: 

AE^ ^ {A\U^^^A)„r, (29) 
where the amplitude U^^'^ is defined by the relation 

(A'lS^^^A) = ~27TiS{EA' - Ea) (A'|C/(2)|A). (30) 

The general proof of Eqs. (l29T l. (|30] | one can find in lfl2ll . 

Insertion of the expressions for the electron (Eq. ( fT9] l) and photon (Eq. (|28] l) propagator into Eq. ( |27| l and integrating over 
time and frequency variables yields 

a4^'-£e(^^^-(-)) ^ (31) 

2tT ^ \ ri2 J AnnA 



InA{ri2)= / 7T-: ^ — ■ — , (32) 

En(l — lO) — Ea + uJ 



where a (i=l,2) are the Dirac matrices acting on the different electron variables. 
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The exact evaluation of the integral ( |32] | in the complex uj plane results in 11231] : 



En\J \ \l3n 



e^l/5.A|n. + [cz(|/3„A|ri2) sin (|/3„^|ri2) - sz (|/3„^|ri2) cos (|/3„A|ri2)] 

\PnA\J \PnA\ 



(33) 



where /3„^ = i?„ — Ea and the notations si(a;) and cA{x) for the integral sine and cosine functions are employed. Then, 
according to Eq. (|26] | the pure one-photon width fJ^^ can be represented in a closed form ll23ll 

2 „ V V V ^12 /a„«A 

where the summation over e and integration over v, remaining in Eq. (fTTT l. is now performed. Note, that the summation in Eq. 
(|34] | extends only over the electron states with energy — m < £"„ < iJ^i. 

The variable /3„A'ri2 is of the order aZ, so that in the nonrelativistic limit /3„a?'12 <C 1. Then expanding sin (|/3„a|'^12) in 
Eq. ( l34l i we can again make contact with the nonrelativistic expression Eq. (fT4l) . or Eq. ( [TtT i for 



VI. DECAY WIDTH VIA THE IMAGINARY PART OF THE LAMB SHIFT: APPLICATION OF THE ADIABATIC THEORY 



In this Section we will apply the Gell-Mann and Low adiabatic formula P 22|] for the energy shift IS.Ea (Lamb shift) of an 
excited atomic state A due to the interaction with the vacuum of the radiation fields 



1. e§-M\S,\A) 
iXbjA = lim —iri—^ — 

-J-oZ {A\S^\A) 



(35) 



for the evaluation of the imaginary part of the Lamb shift. The adiabatic ^-matrix 5,, differs from the usual S'-matrix by the 
presence of the adiabatic (exponential) factor e^ ''!*' in each (interaction) vertex. It refers to the concept of adiabatic switching on 
and off the interaction introduced formally by the replacement Hint{t) — » Hll^^{t) — e^'''*' Hintit). The symmetrized version 
of the adiabatic formula containing S',,(oo, —oo), which is more convenient for the QED calculations, was proposed by Sucher 
ll25ll . The first application of the formula ( [35l l to calculations within bound-state QED was made in ll26ll . In ll26ll it was shown 
how to deal with the adiabatic exponential factor when evaluating the real part of corrections to the energy levels Eq. dSST l (see 
also 1123(1 ). In this paper we will employ the same methods for evaluating the imaginary part of Eq. (l35T l. 

For a free atom (or ion) in the state \A) interacting with the photon vacuum |0-y) (i.e. |A, 0-y) = |A)|0-y) in the absence 
of external fields) the complex energy correction Eq. ( |25] l contains only diagonal S'-matrix elements of even order, since 



(O^jS'^^'' |0-y) = (O^jS'^'''' |0^) = etc. For the separation of the imaginary part of the energy shift AE)^^' of a given order 



:(3) 



2i, it is more convenient to represent Eq. 



in terms of a perturbation series of the form (up to terms e ) [23.1 



AE/ 



lim in 



(36) 



For the adiabatic S,j matrix we used the standard expansion in powers of the interaction constant e 



(37) 



and can separate real and imaginary parts of the matrix elements at any given order of perturbation theory 

{A\Sip\A) = Re{A\Sip\A) +iIm{A\Slp\A) . 
The only one second-order term describes the pure one-photon decay width 



ImAE 



(38) 



(39) 



Arranging aU the terms of fourth order, which describe the pure two-photon decay width including - as we will see below - a 
part of the (one-loop) radiative corrections to the one-photon width, one obtains 



ImAE^f^ — lim i] 

77^0 



(40) 
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where the last two terms result from the expression {A\SrP 

The total width of an excited electron state A (specifying the initial state as \ A, 0^) 
from the imaginary part of the total energy-shift via 



Fa = -2ImAEA , 
respectively, after perturbation expansion of AEa (up to order e^) as 



\A)) should follow (by definition) 

(41) 



= lim —277 



(42) 



The formulas ( [39b , ( l40b and ( l42b will be employed in the next Sections for evaluating the pure one-photon and two-photon 
widths r^'' and F^' , respectively. 



VII. APPLICATION OF THE "OPTICAL THEOREM" 

A. Formulation of the "optical theorem" for the S-matrix elements 

The "optical theorem" is a consequence of the unitarity of the S'-matrix and in the most general case can be formulated as 
follows (see, for example, Il27i1 ). First, we introduce the T-matrix via the definition 

S ^l + if. (43) 

For the diagonal matrix element of Eq. ( |43] ) it follows 

{I\S\I) = l + i{I\f\I) (44) 

and 

Re{I\il ~ S)\I) ^ Im{I\f\I) . (45) 

Here |/) denotes the initial state of the (decaying) system "atom + radiation field". The wave function |/) is assumed to be 
normalized; in our later case of interest it refers to the wave function for the excited atom state and the photon vacuum. As a 
consequence of the unitarity relation S^S = SS^ = 1 for the S'-matrix, the "optical theorem" states 

i (t - f = -f = -f f ^ , (46) 



respectively, for the matrix elements 



2/m(/|f|/)^^|(F|f|/> 



(47) 



The summation in Eq. ( |47] ) runs over the complete set of final states F allowed by the energy conservation law. Formally, 
the state F = / is also included in this summation. The latter circumstance will be important for our further derivations. The 
initial and final states |/) and \F) are by definition eigenstates of the electron and photon number operator. We also mention at 
this point, that the "optical theorem" strictly holds for arbitrary not explicitly time-dependent problems. Specific modifications 
are required in the presence of time-dependent external field, which allow for electron-positron pair creation out of the Dirac 
vacuum. Under such conditions the S-matrix will become nonunitary lEsIl . 

Performing perturbation expansion of the T-(respectively for the S-) matrix one derives immediately from ( l46l ) and (07]) the 
"optical theorem" for the T-matrix elements up to any (even) order 2i of perturbation theory 

F F j<i 

Depending on the physical process (respectively, the scenario) under consideration one has to fix the number of electrons (atomic 
state) and photons (radiation field at zero temperature) in the inital and final state. However, the quantum numbers of electrons 
(A) and photons (k, e) will vary of course. In case of the one-photon decay the summation over F includes the summation over 
the final atomic states A' as well as over the quantum numbers k, e of the emitted photon. In case of the two-photon decay the 
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summation over F includes, apart from the summation over A', also the summation (integration) over the characteristic quantum 
numbers of the two emitted photons. 

Expanding both sides of Eq. ( l45l ) into a perturbation series (with respect to e), we find for arbitrary orders i — 1,2,... the 
relation 



Re{I\S^'^I) = -Im{I\f^'^I) . 
and thus the S'-matrix form of the "optical theorem" corresponding to Eq. ( l48T l 

F F j<i 

Then, collecting the second-order terms in Eqs. ( l47T i and ( fSOl l yields 

-2Re{I\S'^^'>\I) = [(■f'l^^'V) 

F^I 

Again, only nondiagonal matrix elements, like (F|S'(^) |/) contribute. 
Collecting now the fourth-order terms, we find 

- 2i?e(/|^(4)|/> = + |(-^l^^'V)f + H 2Re{I\S'^^^\F)(F\S'^'^^\I). 

F^I Fyil 



(49) 



(50) 



(51) 



(52) 



The last term in Eq. ( |52] | represents, evidently, the radiative corrections to the one-photon width. These corrections were 
evaluated by Barbieri and Sucher via direct evaluation of the corresponding imaginary part of the two-loop Lamb shift (id}. 
Here we will not repeat these calculations within our approach. Note, that the term F = / in the sum over F is absent in this 
contribution, since {I\S'-^Ml) = (^l-S'^^V) = 0- 



B. Application of tlie "optical theorem" to the adiabatic S'-matrix elements 



As indicated above the adiabatic S'-matrix Sri arises after introduction of the adiabatic switching function f{ri) = e^'''*' in 
the QED interaction Hamiltonian. Assuming that no dynamic excitations of the system takes place during switching on and 
off the interaction, the adiabatic S'-matrix remains unitary ll28ll . ||29tl . Moreover, all observables calculated on the basis of the 
adiabatic approach should not depend on the specific form used for the adiabatic factor after the limiting process rj ^ has been 
performed. Therefore, we will apply the "optical theorem" relations ( fSTT i and ( |52] | to the adiabatic formulas ( [39] ), ( l40b and ( l42b . 
In what follows, it will be necessary to fix not only the state of an electron in an atom, but also the number of the photons. 

Then from Eq. ( l39l) . (l42l) and ( |52] | it follows for I = A,0~^ (excited state, no photons) for the pure one-photon width 



.(1) 



lim n 



E 

F=/^A,0.y 



\{F\Si,'^\A,0, 



(53) 



and for the two-photon width 



T^2^^l^^v\^ E \{F\S^'^\A,0^f)f + ^{Re{A,0^\S^'^\A,0.,) 

The remaining term up to order containing radiative-correction effects 

r2'^=\imr, J2 2i?e(v4,0^|S'(i)|F)(F|S'(3)|A,0^) 

will not be considered further, since we are aiming at the two-photon decay. 
Employing now Eq. ( fSTT i, we can rewrite Eq. ( |54] i finally into the form 

r^^' = limJ2 Y \{F\S^'^\A,0,)\\2Y\{A,2,\S^'^\A,0,)\\i ^ \{F'\S(' 

I F^A,0^ 2^ \f'^A,0^ 



(54) 



(55) 



(56) 
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In Eq. ( |56] ) we have to distinguish between the final states (F) and (F') for two-photon and for the one-photon transitions, 



respectively. It is important that the term 



2 



has canceled out in Eq. ( 1541 ). The last but one term in 



corresponding to apparently nonphysical transition A A + 27, but formally present in the sum over F states, will indeed 

cancel out in the final expression (see Section 9). The notation ^ means here the integration over the frequencies of two photons. 

2-, 

In the next Section we will evaluate the one- and two-photon decay widths using Eqs i53[ and ( |56l l. 

VIII. ONE-PHOTON DECAY WIDTH VIA THE "OPTICAL THEOREM" 

We start with the evaluation of the decay width F^' using Eq. ( l53b . First, we evaluate the matrix element {A' , k, e\sij''^ \A, 0-y) 
for the emission of the photon with momentum k and polarization e. This matrix element for the "normal" S'-matrix was 
evaluated in Section 2. The corresponding adiabatic 5^-matrix element reads 

{A',k,e\Sip\A,0^) =e J d^x4>A'{xh^.A;{x)i;A{x)e-''M . (57) 

Now the integration over the time variable yields essentially a representation of the ^-function 

on 

dte^(E.-E,,-u.)t-,M ^ ^ 2nS^{LUAA' - co), (58) 

{uJAA'-^r + V 
— 00 

where lim d„{x) = 6{x). As the next step we perform the integration over the photon frequency. Taking Eq. (|58l l by squre 
modulus, multiplying by ui and integrating, we obtain 

00 

/ r, ^^^^^^^^ - V ( ^ + :^ + ^-rctg ( ^) ] (59) 







Having in mind the limit 77 we can replace Eq. d59] l by 

00 

4^2 / ^ ^ (60) 

It remains to multiply the result by the factor {2n)^^ (see Eq. Q), by the factor {\/2tt) (see Eq. (|2]l) and by the factor 77 from 
Eq. ( l53T l. The matrix element will be the same as in Eq. ^ and we again arrive at Eq. ( fTTT ) with the summation over the electron 
states, lower by energy than the state A: 



A'(E^,<Ea) 



-tkf 



A' A 



(61) 



In the derivation above the manipulations with (5-functions, like in Section 2, have been avoided. Multiplying the result by the 
adiabatic parameter 77 in Eq. ( |53] l plays the same role as dividing the result by the time T in Section 2: the adiabatic factor 77 
has the dimensionality s^^. Note, that in this approach the automatic exclusion (like in Eq. (l34l l) of the transitions to the states 
higher than A in the summation over F in Eq. (l53T l does not occur and we have to refer to the energy conservation law to avoid 
them. 

IX. TWO-PHOTON DECAY WIDTH VIA "OPTICAL THEOREM" IN THE ABSENCE OF CASCADES 

A. Evaluation of the two-photon decay width 

In this Section we will evaluate the two-photon decay width fJ^^ using Eq. ( |56] l. We start with the first term in the curly 
brackets in Eq. ( |56] |. The S'-matrix elements corresponding to the emission of the two photons k, e and fc', e ' look like 

{A',k',e';k,e\Sl,^^\A,0^)a ^ ( d^x^d^x^ {i>A'{xi)^^.,Alf' *{xi)e--'>\''\S{xiX2)l^,,A'j^:*{x2)e~'^\''\il^A{x2))Am 
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where the electron propagator S{xiX2) is defined by Eq. ( fT9] l and the indices a, b correspond to the Feynman graphs Figs 2a, 
2b, respectively. The integration over t2 in Eq. (|62| | results 



2ri 



{uji + Ea — uj) 



(64) 



and similarly looks the integration over ti in Eq. 



2r; 



(65) 



The next step is the integration over uji in Eq. (|62| | which has to be performed in the complex tJi -plane. We have to evaluate the 
integral (considering the iO-perscription in the energy denominator of the electron propagator) 



duji 



[{uji +Ea~ u;)^ + 77'] [(^1 + Ea' + io')^ + [^„(1 - iO) + loi] ' 



(66) 



In what follows in this Section we will restrict ourselves to the nonrelativistic limit of Eqs ( |62] |. (l63T l. since we are most interested 
in the hydrogen case. Then we can fully neglect the sum over the negative-energy states in the electron propagator ( fT9l ). Note, that 
it is possible while we are using "velocity" form for the matrix elements of the photon emission operator In the "length" form 
it would not be the case i24tl . So we can close the integration contour in the lower half-plane, where only two poles are located: 



(jj 



(1) 



Ea + uj + ir] and uj[ = ~Ea' —uj' + irj. In the absence of cascades, the energy denominators (£"„ — Ea + — 
and {En — Ea' + uj' — i'ri)~^ are nonsingular and we can omit the imaginary parts irj in these denominators. Moreover, using 
energy conservation law condition, Ea' + uj' = Ea — uj we can consider both denominators as being equal. Then, collecting 
the factors (l/27r from Eq. ( fT9] l, —2tt from the Cauchy formula for the contour integration in the clockwise direction) yields for 
the amplitude Eq. 



{A\k',e';k,e\S^^^\A,0^)a 



[{ujo -UJ- uj'f + Arf] 
Adding the similar contribution from Eq. ( |63] ) leads to 

(A' 1 5(2) 1^,0^),+, 

r / -1/1 r> 7** l__ 

4?7 



E 



{A'\dAt\n){n\dAl\A) 



En 



Ea+U! 



(67) 



(68) 



[(wo - UJ - uj'Y + W'\ 



WAl^^\n){n\dAt]^A) n) |A) 

^ -Ea+uj ^ ^ 



En - Ea' - UJ 



In the matrix elements of the photon emission operators, unlike the S'-matrix elements, we retain the shorthand notation \A) 
instead of \ A, 0^), since in this case it cannot lead to any misunderstandings. 

Insertion of Eq. ( |68] l into the first term on the right-hand side of Eq. ( |56] l and summation (integration) over quantum numbers 
of the final-state particles yields in the nonrelativistic limit 

r^^^(lst term) = lim 2r] (v^)^^^y E E / J J ^ ^^^^ 



(4r?)2 



[(wo - UJ - uj'Y + 47^2 



{jV\aAt^^,\n){n\aAl_\A) {A'\aAl^\n){n\aAl JA) 

^ En-EA+UJ ^ 



En — Ea' 



The integrations over uj, uj' in Eq. ( |68] | can be performed using exactly the same standard QED procedure as in Section III. We 
first integrated over uj', using the (5-function in Eq. ( |20l i. i.e. over the interval (00, 0), or even over (00, —00) what is actually 
equivalent in this case. The second integration over uj was performed over the interval (0, ujq) (see Eq. ([25])). 
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Let us adopt here this procedure within adiababtic S'^-matrix approach. Integrating Eq. ( [69] ) over uj' according to Eqs 
and ( |60l l leads to 



Lu'duj' 



7r(wo - uj) 



[{uq - CJ - CJ')2 + 4772] 
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(70) 



Then the integration over the emitted photon directions and summation over the polarizations yields in the nonrelativistic limit 
(in the "velocity" form): 



T^^\lst term) = — / uj{ujq ~ uj)duj ^ \iUik{uj)) j^, 

n i,k—l 



(71) 



I A' A 



^ {A'\p,\n){n\pk\A) ^ {A'\pk\n){n\p,\A) 

^ En-EA+LU ^ En-EA'-OJ ' 



(72) 



where pi = {p)i. For A = 2s, A' = Is Eqs. (TtTI i. ( |72] | go over to Eqs (I22b-(l24li. if we use again the quantum mechanical relation 

B. Cancellation of singularities 

Apart from the first term in Eq. ( |56] l there are two additional terms which contain the singularities with respect to the adiabatic 
parameter 77 in the limit 77 ^ 0. In this Subsection we will show that these singularities exactly cancel. We start with the last 
term in the right-hand side of Eq. ( |63] |. In the absence of cascades, i.e. when there are no energy levels between the initial 
state A and the final state A' , the only term in the sum over F' \s F' — A' ,1^. Then, after summation over the emitted photon 
polarization and integration over th emitted photon directions, repeating the derivations performed in Section 8, we obtain the 
result 



\F'^A,() 



{F'\S'^^^\A,Q^ 



iimi(r«)^ 



(73) 



This divergent term can be canceled only by the "unphysical" contribution 2 ^ {A, 2~f\S)^ \A, 0-y) in the right-hand side of 

27 

Eq. (|56] |. The latter one looks exactly like Eq. ( |69l ) if A' is replaced by A. Setting n — A' in the sum over n in the expression 
2(A, 2-y|S'^^^ I A, O7) will give the same set of the matrix elements as in Eq. ( |73] ). This contribution also appears to be divergent 
like rj^^ in the limit 77 — > and cancels the divergency Eq. (|73] |. 
To demonstrate this, we write down the expression 



lim 277 y (A,2^|S'(2)|A,0^ 

27 



^lim 277 (27r)2— j dVdV j u^diu J uj'du;' 



(477)2 



[{luo - uj - uj'y + 4:rf 



{A\dAl JA'){A'\dAt\A) {A\dAt\A'){A'\aAt JA) 



Ea' — Ea+ uj + if] 



Ea' — Ea — uj + ii] 



(74) 



In this case, unlike (I68I I, we now keep i77 in the energy denominators, in order to keep trace of all the divergences. The integration 
over UJ is exactly performed like in Eq. (fTOl l with the result 

lim 277 Y: |(A2,|^f |AO,)f = lim 2VT^,T.J ^^^^-'4 | (A|dly f | 1^) 



UJ duj 



1 



1 



Ea' — Ea + uj + ii] Ea' — Ea - uj + irj 



(75) 
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After summation over the polarizations and integration over the emission angels and transforming the expression in 



get 



lim 2ri 



E 

27 



(A 2^1^(2) 1^,07) 



= — lim — 



(76) 



where the function 



F{V) 



(11) 



remains to be calculated. In this case we have to evaluate the integral over lo in the same way as for deriving the expression ( |73] |. 
i.e. integrating over the frequency interval (0, 00). However, it will be more convenient to extend formally the integration over 
the entire interval (cx3, —00) as it was done, for example, in [29]. Then the integration can be performed in the complex plane: 



to dio 



(78) 



The denominator in Eq. (fTST l contains 4 single poles: one pole in each quadrant of the complex plane. However, we have to 
recall that two of these poles originate from the negative (i.e. unphysical) oj values. Therefore we have to omit their contribution. 
Then the evaluation of the integral ( fTST i yields with the aid of Cauchy's theorem 



2lTl TT 

Airj 2ri 



(79) 



Inserting this result into Eq. ( |76] | we find that this result exactly cancels the divergent contribution Eq. (l73T l. What concerns 

" (2) 

the terms with n A' in the sum over n in the expression for 2{A, 2j\S,j \A, 0^), we will avoid them referring to the energy 
conservation law, exactly like we did it in the Section 8 for the one-photon transition. Aiming to cancel the contributions from 
Eq. (|73]l and from 2{A, 2^\S^^^^\A, 0^) we have to treat both contributions in the same way. 

in the absence of cascades cancel each other and the two-photon width is 



Thus, the second and the third terms in Eq. 
given exclusively by the first term in Eq. (|56] |: r^'(lsi term) — F^'. The expression ( f72] i for F^' coincides exactly with the 
standard QED expression for the two-photon decay width in the absence of the cascades. In case of the 2s-level in hydrogen we 
return to the expression (|22]) for n — 2. 



X. TWO-PHOTON DECAY WIDTH VIA "OPTICAL THEOREM" IN THE PRESENCE OF CASCADES 

The derivation of the expression for the two-photon decay width in the presence of the cascades does not differ, in principle, 
from the derivation of Section IX. The expression ( |69] | holds in this case as well. We will assume that there exists only one 
cascade channel (for example, 3s — 2p — Is in case of the decay of 3s-level in hydrogen). The only difference is the existence 
of an additional resonance by n = i? in the sum over n in Eq. ( |69] l: uj = Ea — En {R is the resonance state). Due to the 
energy conservation this implies also the existence of another resonance (lower branch of the cascade): uj' = E^ — Ea'- Now 
Eq. ( |73] ) contains two divergent (like rj^^ hy rj ^ 0) terms: by F' = A', 1^ (this divergency is the same as in case of the 

absence of cascades) and by F' = R, 1^ (this is the additiona divergency connected with the existence of the cascade). The 

^ (2) 

former divergency is compensated by the term n = A' in the "unphysical contribution" 2 (^2^ IS",, \A, Oj) and the latter one is 
compensated in the same way by the term n = i? in the expression for 2(^2^15(2)^^1^, 0-,,). However, in the presence of the 
cascade a third divergency arises directly in the summation over n (for n = R) in the expression (f74l) . This divergency cannot 
be canceled by any counterterm since it is proportional to a special product of matrix element: |(yl'|pi|i? >< No 
counterterm in Eq. ( |56] | contains such a product. This remaining divergency can only be removed by taking into account (at 
given order of perturbation theory) the radiative corrections to the level width. Equivalently, one may introduce the level width as 
it was done in Section IV within the framework of standard QED. Thus, evaluating the two-photon decay width in the presence 
of the cascades, we again return back to the same expressions and the same problems which were discussed in Section IV. 

Our analysis of the evaluation of the two-photon decay width in the presence of cascades via the imaginary part of the Lamb 
shift, performed in Sections IX, X contradicts to Jentschura's "alternative" approach iflTll - lfTgll . First, the expression ( fTTT i contains 
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the square modulus | (i7ife(a'))^,^|^, not the squared matrix element as used in Jentschura's calculation. Second, 

there is no chance to compensate the cascade divergency in Eq. ( |69] l within the evaluation of the imaginary part of the Lamb 
shift via the adiabatic S'-matrix approach. In Jentschura's evaluation of the two-photon decay width in the presence of cascades 
the integral remains finite without introduction of the level width. We claim that the cascade contribution to the two photon 
decay width remains infinite without an introduction of the level widths in the energy denominators (via partial resummation of 
radiative corrections). 



XI. EVALUATION OF THE TWO-PHOTON WIDTH FOR THE 3s LEVEL IN HYDROGEN 

As a consequense of the studies presented in Sections 5-10, we return now to the standard QED expression for the two-photon 
decay of 3s-level in hydrogen (Eqs (I22li-(|24]| with n = 3) and to the prescription for the employment of the formulas (l22Tl-(l24ll 
given in Section 4. Inserting Eq. (|22] | into Eq. (|25] | and retaining only the resonant term within the second and fourth frequency 
intervals, will yield the cascade contribution to the total two-photon decay rate of the 3s-level. Taking the ratio to the total 

width of the 3s-level T^s we wil obtain the absolute probability or branching ratio M^g^."^/'"''^' /^Ss = b^^s-2p^is f^i" cascade 

transition. The contributions to fog^'^f.^p-'is from the intervals (I), (III), (V) are assumed to be zero. The cascade contribution of 
the 3s-level results (in the "length" form) 



3s;ls 



277r 



(11) 



{R3s{r)\r\R2p{r)){R2p{r')\T'\RUr')) 



— + — 



dhJ 



(80) 



(IV) 



{RUr)\r\R2,{r)){R2,{r')\r'\RUr')) 



— Els ^ ijj — f 



dio. 



According to the discussion in Section 4 the "pure" two-photon decay probabilities within each interval, defined in Section 4, 
look like 



4 
27^ 



1 2 



dtj, cj G II 



277r 



1 2 



dw, w e IV 



277r 



W^(W0 - CJ)^ [S'ls;3s(^) + 'S'ls;3s('^0 - W)] ^tj, W G I, III, V . 



(81) 
(82) 
(83) 



HereS'j^.f3^^(w) is the expression with the n — 2 term being excluded. 

Unlike cascade, all the intervals contribute to the "pure" two-photon transition. The branching ratio for this transition 3s 
27 + Is appears to be 



(84) 



It remains to introduce the interference contribution. This contribution we consider only for the 2nd and 4th intervals. The 
corresponding frequency distribution functions are given by 



(inter) 1 _ 4^3(^0 - Lof 



277r 



Re 



(i?3s(r)|r|j?2p(2r))(fl2p(rO|r-|Ei3(/)) 
E2p — E'is + — |r 



duj (85) 



277r 

and branching ratio results as 



{Rzs{r)\r\R2p{2r)){R2p(r%r'\Ris{r')) 

E2p — Els ~ ~ f r 



5i^;3s(w) + S'|^f3^,(wo-w) 



duj (86) 



, (inter) 
"3s;ls 



2r 



3s 



2r 



3s 



(") 



(IV) 



(87) 
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The results of our calculations are presented in Table 1. It is convenient to define the size Aw of the second interval as 
multiples I of the widths F, i.e. Aoj — 2ir and for the fourth interval as Aw — 2ff2p, respectively. In Table 1 numbers are 
given for different values of / ranging from I ~ 10^ up to Z ~ 10^. The upper bound of interval II equals lui + IT = + W (in 
a.u.), while the lower bound of interval IV equals uj2 — 2p = | — l^2p- The different lines of the Table 1 present branching 
ratios and transition rates of the "pure" two-photon and "interference" channel, respectively. For the more detailed analysis the 
contributions of the "pure" two-photon transition rate for the each frequency interval are also compiled. The branching ratio 

and the transition rate for the cascade contribution can be obtained from the relation b'"^^22^°is + ^2s"ir^^'' + ^3s"is ' = 1- This 
relation is sutisfied with high accuracy since the only decay channel neglected is the very weak direct 1 -photon Ml transition 
3s Is + 7. From the Table 1 we can draw the following conlusions: as in the case of the HCI ifTTIl . the "pure" two-photon and 
cascade contributions to the total decay rate appear to be inseparable. Changing the interval size Aw, we obtain quite different 

values for dw!^^^''^''^ ranging from 202.16 s^^ (for I = lO'^) up to 7.9385 s^^ (for / = 1.00256 • 10^). 

Moreover, in our calculations - depending on the size of the interval - the interference contribution also can become quite large, 
comparable in magnitude with the "pure" two-photon contribution. Thus, we demonstrated that even the order of magnitude of 
the "pure" two-photon decay rate for the 3s-state in hydrogen can not be predicted reliably. 

Earlier the result 8.2196 s^^ for the "pure" two-photon decay of the 3s-level was reported in [jst] and confirmed in i^. 
However, as it was pointed out in |14] in both papers f^, ^ the summation over the intermediate states was not performed 
properly. The "nonresonant" contribution 10.556 s^^ deduced in 114.1 . which plays the role of the "pure" two-photon decay rate 
is well within the range of our values given Table 1. However, the result 2.08 s^^ obtained for the "pure" two-photon decay rate 
in lfT6ll is in strong contradiction with the present analysis. 

Very recently, a paper |30] did arrive where both the standard QED approach, based on the line profile theory (||9|]-|fTT|]) and 
the "alternative" approach based on the two-loop Lamb shift theory ( UtII - IIIQII ) were applied to the calculation of the two-photon 
transition in hydrogen. A reasonable agreement between the two methods was found. However, from the derivations in our 
present paper it follows that the employment of the Lamb shift imaginary part gives exactly the same results as the standard 
QED approach. The difference between the "standard" and the "alternative" methods is due to the use of the squared amplitude 
instead of square modulus in Jentschura's calculation. To our mind this replacement is unacceptable and cannot be justified 
within QED. 



XII. CONCLUSION 



In this paper we developed a method for the calculation of the two-photon decay rates, based on the evaluation of the imaginary 
part of the Lamb shift with employment of the adiabatic 5-matrix theory and the "optical theorem". We have shown that the 
results of such calculations coincide exactly with the standard QED approach also in the presence of cascades. 

We demonstrated that a strict separation of the "pure" two-photon and cascade contributions for 3s-level decay in hydrogen 
is impossible. Moreover, we show that even the approximate separation of these two decay channels cannot be achieved with an 
accuracy, required in modern astrophysical investigations (i.e. at 1 % level) of the recombination history of hydrogen in the early 
Universe. 

As a possible solution of the problem with respect to astrphysical needs, we would suggest to rewrite the basic evolution 
equation for the number of the hydrogen atoms in a certain excited state (i.e. Eq. (2) in 1.15.1 ) in a way which does not distinguish 
"pure" two-photon decays and cascades. 
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A' 



Fig. 1. The Feynman graph corresponding to the photon emission in an one-electron atom. The double solid line describes the electron 
in the firld of nucleus (Furry picture), the wavy line with the arro at the end describes the emitted photon. The indices A and A' refer to the 
quantum numbers of the initial and final states of an atom, cj denotes the frequency of the emitted photon. 
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Fig. 2. The Feynman graphs corresponding to the two-photon transition A ^ A' + 27. By lu, uj' we denote the frequencies of the emitted 
photons. 
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Fig. 3. The frequency distribution d,W.^^^-{^ for the total two-photon transition 3s ^ Is + 27 including cascade and "pure" two-photon 
transitions as functions of the frequency (in a.u.). The transition rate divided by (a is the fine structure constant) is plotted versus the 
frequency within the interval [0, ujo]. The boundaries for the frequency intervals I-V are also indicated as vertical lines. 




A A 



Fig. 4. The Feynman graph corresponding to the one-loop electron self-energy. The imaginary part of the contribution of this graph to the 
energy level shift AEa represents the one-photon width Fa of the level A. The wavy line denotes here the virtual photon. 



TABLE I: Branching ratios and transition rates (in s ^ ) for the different decay channels for the decay probability of the 3s level with different 
frequency interval size (I). 



I 


10^ 


10^ 


2.5 ■ lO-""' 


5 • lO '' 


10" 


1.5 ■ 10** 


4.53 ■ lO'' 


1.00256 ■ 10^ 


^(pure27) 


3.2003 • 10-=^ 


3.5091 • 10-'^ 


1.6270 • lO"'' 


1.0239 ■ lO"'' 


7.6765 • 10^^ 


7.2201 ■ lO"'^ 


9.1487 • lO^'^ 


1.2567 • lO"'' 


|y(pure27) 


53.054 


7.0547 


3.5743 


2.1898 


1.27737 


0.85130 


2.4979 • 10"^ 







0.006247 


0.06247 


0.15614 


0.31201 


0.62183 


0.92718 


2.4666 


3.9810 


TI/{purc27) 
VV jlj 


95.536 


7.8778 


2.7928 


1.4517 


1.0457 


1.0031 


0.86005 





TJ/(puro27) 


0.006185 


0.061847 


0.15458 


0.30890 


0.61569 


0.91813 


2.4523 


3.9575 




53.561 


7.1101 


3.5999 


2.2056 


1.2886 


0.861254 


3.1665 • 10"* 





^/{puro27) 


202.16 


22.167 


10.278 


6.4680 


4.8492 


4.5609 


5.7792 


7.9385 


^(inter) 


-1.4342 • 10"'' 


-1.4343 • 10"* 


-3.5852 • 10"* 


-7.1665 • 10"** 


-1.4302 • 10"^ 


-2.1376 • 10"^ 


-6.0829 • 10"^ 


-1.0459 • 10"'^ 


pj^ (inter) 


-0.0090599 


-0.090602 


-0.22647 


-0.45270 


-0.90346 


-1.3503 


-3.8426 


-6.6067 



